Roton confinement in trapped dipolar Bose-Einstein condensates 



M. Jona-Lasinio, K. Lakomy and L. Santos 
Institut fur Theoretische Physik, Leibniz Universitdt, 30167 Hannover, Germany 

(Dated: April 5, 2013) 

Roton excitations constitute a key feature of dipolar gases, connecting these gases with superfluid 
helium. We show that the density dependence of the roton minimum results in a spatial roton 
confinement, particularly relevant in pancake dipolar condensates with large aspect ratios. We show 
that roton confinement plays a crucial role in the dynamics after roton instability, and that arresting 
the instability may create a trapped roton gas revealed by confined density modulations. We discuss 
the local susceptibility against density perturbations, which we illustrate for the case of vortices. 
Roton confinement is expected to play a key role in experiments. 
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Dipolar gases have attracted a growing attention in 
recent years. Quantum degenerate gases of magnetic 
atoms as chromium pQ, dysprosium [2] and erbium [3] 
have already been realized. Moreover, the preparation of 
heteronuclear molecules in their ro-vibrational ground- 
state [4] opens the path for the creation of a degenerate 
gas of polar molecules, a goal currently pursued by var- 
ious groups worldwide [5 . Rydberg atoms provide yet 
another possible realization of a highly polar gas [6]. 

The rich physics of dipolar gases arises from dipole- 
dipole interactions (DDI) [7 . Dipolar Bose-Einstein con- 
densates (dBECs) feature a geometry-dependent stabil- 
ity [8] and a peculiar dispersion of the elementary exci- 
tations. Nonpolar BECs present the usual Bogoliubov 
spectrum, with a linear (phonon) dispersion at low mo- 
menta, and a quadratic dispersion at large momenta [9]. 
Under proper conditions, dBECs present a dispersion 
minimum at intermediate momenta [lOj [11] resembling 
the roton minimum of superfluid He [T2] . 

Roton excitations are crucial in He, reducing the criti- 
cal superfluid velocity [13^ and leading to density modula- 
tions at defects [T5HT9] . Remarkably, similar effects have 
also been predicted in dBECs [T41I20H22] . Moreover, the 
roton minimum is crucial for the stability of a dBEC. 
When the minimum reaches zero energy, the dBEC be- 
comes unstable against finite- momentum excitations (ro- 
ton instability) [TQl [TTJ [23] , which differs fundamentally 
from the usual phonon instability. 

Roton properties in He may be controlled by means of 
pressure [24 . Similarly, the roton minimum in dBECs 
depends on interactions and therefore on density. In this 
Letter we demonstrate that this density dependence leads 
to a spatial roton confinement in trapped dBECs. Roton 
confinement has been hinted in recent numerical calcu- 
lations p~4j [25] and resembles that of rotons at vortex 
lines in He [26]. We show that roton confinement is cru- 
cial to understand the roton instability in pancake traps. 
We also discuss how arresting the instability may create 
a trapped roton gas. We finally analyze other conse- 
quences of the roton confinement, such as local suscepti- 
bility against density perturbations. 




FIG. 1: (Color online) (Top) Local spectrum of a BEC of 
2 x 10 5 Er atoms uo z = 2tt x 1 kHz and A = 40. Note a 
minimum (dark red region in the top figure) in both space 
and momentum. (Bottom) Change of the roton minimum for 
two different radial positions. 



We consider a dBEC of N bosons of mass m and (elec- 
tric or magnetic) dipole moment d oriented along z. The 
dBEC is in a pancake harmonic trap V t (r) of frequencies 
uj in the xy plane and uo z = Xu along z, with A ^> 1. The 
dBEC wave function 0(r,£) obeys the nonlocal Gross- 
Pitaevskii equation (GPE) [7], 
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where g = Aith 2 aN/m characterizes the short-range in- 
teractions, a is the s-wave scattering length, Vdd(r) = 
^pr(l - 3 cos 2 8) is the DDI potential, 6 is the angle be- 
tween r and the z axis, and J <i 3 r|(/>(r, t)\ 2 = 1. 

We first consider uj = [27], briefly summariz- 
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ing the results of Ref. [TP] . The ground-state wave- 
function is <fro(z) exp (—ifjit/h), where /i is the chemi- 
cal potential and 4>o(z) fulfills a ID local GPE with 
a regularized coupling constant g + g^, with = 
87rNcP/3. Assuming a transverse Thomas- Fermi (TF) 
profile (j>o(z) = \/uq (1 — z 2 /L 2 ), with no the peak 
density, one obtains \i — (g + ga)^. Excitations 
of energy e are evaluated by substituting (j)(p,z,t) = 
e -ifit/h |0 o ^) + z )e- iet /^ - v (p, z)*e ie ^] (with 
p = (x,y)) into Eq. ([!]), and keeping only linear terms 
in u and v. The excitations have a defined in-plane mo- 
mentum q, and hence f±(p,z) = [u(p, z) ±v(p, z)] = 
f±(z)e LCi ' p . The lowest eigenenergy obtained from the 
Bogoliubov-de Gennes (BdG) equations for each q builds 
the dispersion e(q). Interestingly, e(q) may present a 
roton-like minimum at intermediate q. Assuming qL ^> 1 
and nE(q)\P - 2|/(/3 + 1) < H 2 uj 2 z (with /3 = g d /g and 
i£(g) = h 2 q 2 /2m) we obtain an approximate expression 
of the dispersion, e\(q, p) E(q) 2 - G(f3)E(q)p + ft 2 u; 2 , 
with G(/3) = 3p^y^^j7 as we U as the associated eigen- 
state, f+(z) c± 0o (z). The expression e^(g, /i) agrees well 
with the numerical e(g) and shows that the roton depth 
depends explicitly on \i and hence on density. 

Interesting insights about the roton physics for uo > 
are provided by the concept of local spectrum e(g, p), 
which we introduce here. We compute from Eq. ([!]) the 
ground-state no(r) = |</>o(r)| 2 , which in the TF regime is 
approximated by no(r) = ho (l — p 2 /R 2 — z 2 /L 2 ). We 
obtain for each p the z-profile tiq D (z) = no(r)/ J dzno(r), 
and the local chemical potential pi(p). Solving the corre- 
sponding ID BdG equations [10] we obtain €h(q, pi(p)), 
which approximates the local spectrum e(q, p) consis- 
tently with the LDA. The local chemical potential pi(p) 
decreases quadratically with p, and hence at the trap 
center (p = 0) the roton energy is lowest: e r /hu z = 
e(q = q r ,p = 0)/hw z ~ ^1 - (G(/3)^(0)/2^) 2 . The 
local spectrum presents a minimum both in momentum, 
at q r l z ~ ^G(f3)pi(0)/huj z (note that q r l z ~ 1, with 
l z = ^h/muj z , see [10]), and in space, at p = (Fig. [lj). 
For larger p the minimum becomes shallower and even- 
tually disappears. Hence, remarkably, the inhomoge- 
nous BEC density results in a spatial roton confinement. 
Around the minimum: 
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where the effective roton mass 
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and the effective harmonic frequency uj* = cj,_ 
define the roton localization length I* = ^/h/m*uu* = 
2 1 / 4 (i?/g r ) 1 / 2 . Note that - and hence Z* <C 

i? for A > 1. Moreover, q r U ~ a/#7^ > 1 if A > 1, 
which justifies the use of LDA above. 

We now calculate the localized roton wave functions 
using the LDA and the formalism of Ref. [10] . Assuming 



Z* <C i?, we approximate f+(p,z) ~ F(p)0o(p, z), where 
F(p) has a narrow momentum distribution F(q) centered 
around q r with a width Sq oc 1//* <C g r . The finite width 
of -F(q) must be now considered and as a result F(q) ful- 
fills the eigenvalue equation e 2 F(q) ~ [e r +i^] 2 F(q), with 
H = 2^r(^ ~ Qr) 2 — ^^*^*Vq. Interestingly, i7 resem- 
bles the Hamiltonian of a trapped BEC in the presence 
of spin-orbit coupling [29 , 30], where the Rashba-like dis- 
persion ~ (q — q r ) 2 acts as a ring-like potential in q-space. 
The eigenfunctions of i^, T] niS (q)e %Sip / \fq^ fulfill 
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For g r Z* ^> 1, we expand around q ~ q r 
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obtaining the 



eigenenergies E ni3 /huj* ~ (s 2 - l/4)/(2(<? r /*) 2 ) + n + 
1/2, characterized by the angular momentum s around 
the Rashba-like ring, and the radial harmonic excita- 
tions, n, with frequency uj*. The lowest roton states 
have n = 0, being in real space of the form ip s (p) ~ 
e isif e - P /' 2 Lj s (q r p^ w ith J s the Bessel function. 

The localized states ijj s are crucial for the dynam- 
ics following roton instability. We consider a stable 
BEC prepared with an initial scattering length ai > a c . 
The critical a c for the onset of instability depends non- 
trivially on A and g& [7 ; for the parameters of Fig. [2] it 
is a c = 8.48 ao, with ao the Bohr radius. We are inter- 
ested in the instability dynamics after a sudden quench 
to af < a c . Unstable modes lead to a modulational in- 
stability characterized by a growing density perturbation 
5n(r,t) = n(r, t) — no(r) oc ^no(r)3l(f_(r)). As for the 
uj = case, /_ ~ (E(q r )/e)f+ ex ^no(r)F(p). Hence 
the most unstable roton modes result in a localized mod- 
ulation at the trap center, Sn(r) /rio(r) oc 3?(^ s (p)). 

The most unstable mode is ^o(p) 7 but other modes 
may contribute to the instability due to the small en- 
ergy difference between low-s levels (~ huj*/(q r l*) 2 ). As 
a non-trivial consequence of that, the density pattern 
that develops after the quench is influenced by the ini- 
tial population of the excitations exponentially ampli- 
fied during the destabilization. We mimick this depen- 
dence by considering a small initial seeding 0(r, t = 0) = 
(/>o(r)e* % ( r ), where 0o is the ground-state calculated for 
a^, and |x(r) |/7r < £ is a random phase sampled from a 
homogeneous uniform distribution with a variable ampli- 
tude £. Although this allows us to discuss the possible 
collapse scenarios, the actual amplitude of the initial fluc- 
tuations depends on a$ and on the temperature T, and 
its analysis lies beyond the scope of this Letter [33 . 

Figures |2] show the results of our simulations of Eq. ([!]) 
for an erbium dBEC [34]. For a small initial population 
of the unstable modes, the modulation instability pro- 
ceeds at a sufficiently slow pace and ipo(p) dominates. 
As a result, a localized pattern of concentric rings devel- 
ops (Fig. [2|i), in excellent agreement with 5n(r)/no(r) oc 
^o(p), see Fig. EH. The corresponding momentum distri- 
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FIG. 2: (Color online) Roton instability for 10 5 Er atoms, with uo z = 2n x 450 Hz (l z ~ 0.3 /mi), A = 30, en = 8.49a > 
a c — 8.48ao and a/ = (see text), (a) Column density n 2D (p) = J n(r)dz showing concentric rings (s = 0) formed after 
t = 19 ms for small initial fluctuations (£ = 10 -10 , see text), (e) Modulational instability after t = 15.5 ms consisting 
of several s states for large initial fluctuations (£ = 10 -6 , see text), (b) and (f): momentum distribution of (a) and (e) 
respectively (we have suppressed the large peak at q = y/q% + q$ = 0). (c) and (g): post-collapse dynamics after t = 23 
ms for (a) and t — 19.5 ms for (e), respectively, (d) Radial cut of (a) (green crosses) and theoretical column density of 
the form n 2D (p) oc (1 - p 2 /R 2 ) 3/2 (l + AJ (q r p)e- p2 /21 * ) (solid black line), with R ~ 51 Z z , A ~ -0.4, ^ r / z - 1.25, and 
/* = 2 1 / 4 (R/q r ) 1 / 2 ~ 7.6l z . The value of g r is calculated using the local spectrum picture, (h) Remnant atoms for the case (a) 
(solid red line) and (e) (dashed blue line). 




bution presents a Rashba-like ring (Fig.J^j}) [35]. In con- 
trast, for larger initial fluctuations the pattern growth is 
too fast to select only, and the created density pattern 
results from a (shot-to-shot dependent) linear combina- 
tion of different -0 S , being characterized by a superpo- 
sition of eccentric collapse centers (Fig. [2^). The corre- 
sponding momentum distribution still presents a ring-like 
structure but with an azimuthal modulation arising from 
the linear combination of various ^ s (Fig. [2^). 

The global (phonon-like) collapse studied in chromium 
and erbium dBECs [3j [36] results in large atom losses 
and in a <i-wave pattern in time-of- flight (TOF) pictures. 
Remarkably, the roton instability discussed above leads 
to a very different collapse dynamics. Three-body losses 
become crucial in the collapse dynamics and we included 
them by adding -ih^fN 2 \(j)(r, t)| 4 0(r, t) [36] to Eq. 0, 
with a loss rate L3 = 10 _28 cm _6 s _1 . The concentric 
rings appearing in Fig.[2^i eventually undergo a sequential 
collapse and azimuthal instability, starting from the in- 
ner (denser) ones towards the outer ones (Fig. [2^) leading 
to a step-like atom number decrease (Fig. [2J1) . The su- 
perimposed eccentric collapse centers appearing in Fig. [2^ 
lead to a complex post-collapse behavior with character- 
istic mutually interfering jets expelled out of each local 
collapse center (Fig. [2^). In this case the atom decrease 
is smooth (Fig.[2]i). In both scenarios TOF pictures lack 
the (i-wave symmetry, presenting, due to the local char- 
acter of the collapse, a large peak at q = 0, with an addi- 



FIG. 3: (Color online) Same parameters as in Fig. [2p but 
arresting the collapse after 16 ms with a quench of a up to 
8.55 ao > a c . Snapshots of the central region right after the 
quench (a) and 80 ms later (b). Note a clear density modu- 
lation confined at the trap center. 



tional ring feature during the first collapse stages (see a 
thorough TOF analysis in the supplementary material). 
Moreover, the local collapse results in small atom losses, 
compared to the large losses in global collapses. Roton 
local collapse can be hence distinguished from the phonon 
global collapse by the atom losses and the TOF images, 
even without accessing in-situ dynamics [37] . 

Interestingly, the roton instability may be used to cre- 
ate a confined roton gas. The density pattern obtained 
after quenching down to df < a c may be interpreted as 
the growth of the roton population. At the initial stage 
of the instability, the roton population is small compared 
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FIG. 4: (Color online) Vortex lattice for a BEC of N = 10 5 
Er atoms at the threshold of the roton instability. Parameters 
as in Fig. [I] with a rotational frequency 0.3 u. 



to N and we may neglect condensate depletion or roton- 
roton interactions. Once the roton gas is populated we 
stabilize it by quenching up to a > a c . Since m* and uj* 
do not vary significantly around the instability threshold, 
the created rotons (density patterns) remain confined at 
the trap center after the quench, as shown in Fig. [3] 

A deep roton minimum induces close to a perturba- 
tion a large susceptibility against the formation of density 
modulations with the roton wavelength. This well-known 
effect in He [IMS] is also relevant for dBECs [21J. The 
local spectrum picture implies a spatially dependent sus- 
ceptibility, which we illustrate for the case of vortices [38] . 
Vortex cores present a crater-like shape in the case of a 
deep roton minimum [20j [21] , that disappears when the 
minimum is shallow or absent. Therefore, vortices at dif- 
ferent positions in a trapped BEC present a different core 
profile. This is shown for a vortex lattice in Fig. [4] where 
we depict the ground-state of an erbium BEC rotating 
around the z axis with an angular frequency Q = 0.3 uo. 
Note that core modulations at the trap center disappear 
for vortices at the boundary. 

In summary, we have shown that an inhomogeneous 
trapping in pancake dBECs with large aspect ratios leads 
to a spatial roton confinement which is crucial to under- 
stand the roton instability. The roton dispersion has not 
yet been observed experimentally, being currently a ma- 
jor goal pursued by several groups. Roton confinement 
is expected to play a key role in these experiments, since 
harmonic traps are typically employed and large aspect 
ratios are needed to study the roton dispersion. In addi- 
tion to the local susceptibility discussed above [38] , roton 
confinement should be carefully considered when mea- 
suring the critical superfluid velocity, performing Bragg 
scattering [39], or analyzing finite temperature physics, 
which may be very interesting since the thermal roton 
cloud is expected to localize at the center of the trap. 
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dation, and the DFG (SA1031/6 and Excellenzcluster 
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SUPPLEMENTARY MATERIAL FOR "ROTON 
CONFINEMENT IN TRAPPED DIPOLAR 
BOSE-EINSTEIN CONDENSATES" 

DERIVATION OF THE LOCALIZED ROTON 
WAVE FUNCTIONS 

Homogeneous xy case (cj = 0) 

We briefly summarize in this subsection the main for- 
malism developed in Ref. [10 j. For lj = 0, the ground- 



state is of the form <fio(z)e % ^ l l h ^ where \i is the chemical 
potential and 



2 m 



m 



V 2 + -uj 2 z z 2 + (g + g d )\M*)\ ~ » } M*) = 0, 

(4) 

We assume (g + g^) > 0. For \i ^> hu z the condensate 
presents a Thomas- Fermi (TF) density profile, |(/>o(z)| 2 = 
no(l — z 2 /L 2 ), with no = /i/(g + go) being the maximum 
density, and L = ^/2ji/muj 2 the TF radius. 

Excitations of energy e are eval- 
uated by substituting 0(p, z, t) = 
e -ifit/h [0 o( ^ + u ^ z ^ e -iet/h _ v ^ z y e iet/h] ( with 

p = (x,y)) into the 3D non-local GPE (Eq.(l) of our 
Letter), and keeping only linear terms in the Bogoli- 
ubov amplitudes u and v. For uj = 0, the excitations 
have a defined in-plane momentum q, and we may 
write f±(p,z) = [u(p,z)±v(p,z)] = f±(z)e i *>>. The 
Bogoliubov-de Gennes equations (BdG) read: 



2m 



' dz 2 



e/+ = H k inf- + H int [f-], 



f+ = H k i n f '+, (5) 
(6) 



with 

H int [f_] = 2(g d +g)f_(z)\Mz)\ 2 ~ ^/2)g d qcf> {z)x 

rOO 

x / dz' f.(z')Mz')exp(-q\z-z'\). (7) 



For each q we calculate the eigenenergies. The lowest 
one, e(q) builds the dispersion law. 

The most interesting behavior occurs for qL ^> 1. We 
may now introduce f+(z) = W(z)<fro(z). Expressing /_ 
through W from Eq.([5|, we substitute it into Eq.Q and 
integrate over dz' in H int [f_] as the main contribution 
to the integral comes from a narrow range of distances 
\z' — z\ ~ l/q. This yields (here x = Z /L) 
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X X >d X 2 



2(l + (3)J /v d X 



x- 



d 



W( X )- 



E(qf 



2/3-1 
1 + 13 



»E(q)(l- X 2 ) 



3H 2 uj 2 
2(1 + 0) 



W( X ) = 0, (8) 



where E(q) = h 2 q 2 /2m. Here we omitted terms of the 
order of E(q)h 2 uj 2 / ji and ft 4 cj 4 //i 2 , since they are small 
compared to either h 2 uj 2 or E(q) 2 . For each mode of the 
confined motion (along the z direction), the solution of 
Eq.([8| can be written as an expansion series in Gegen- 
bauer polynomials C*(x), where A = (4+/3)/2(l+/3), and 
n > is an integer. For fiE(q)\2(3- l|/(l + /3) < H 2 uj 2 the 
lowest solution is of the form W(x) — l + Xln>o a nCn(x)? 
with small amplitudes a n of Gegenbauer polynomials of 
higher order, while the lowest-eigenenergy (and hence the 
dispersion) is given by 



e 2 (q^)=E(q) 2 -G(f3)E(q)p + h 2 



(9) 



with G(/3) 



QS-2)(5/3+2) 
3(l+/3)(2/3+l)' 



Trapped xy case 



We consider in the following uj > 0, with A = uj z /uj ^> 
1. Due to the x?/ confinement, we cannot employ any- 



more a plane wave basis for the Bogoliubov amplitudes, 
as we did above. However, we assume that: (i) the Bo- 
goliubov amplitudes are, in momentum space, strongly 
peaked around a momentum q ^> l/R; and (ii) the Bo- 
goliubov amplitudes are localized at the center of the 
trap, in a region with a radius much smaller than R. 
These assumptions must be checked self-consistently, but 
we anticipate that they are fulfilled for a sufficiently large 
aspect ratio A ^> 1. If this is the case, we may proceed as 
for the homogeneous case, but taking into account that 
(a) the chemical potential is now p dependent, and (b) 
the Bogoliubov amplitudes are not plane waves anymore. 
Assuming the above conditions are met, we may write 
f+(p,z)czF(p)W(x)MP,z) (now x 2 =p 2 /R 2 +z 2 /L 2 ), 
where the function F(p) has a narrow momentum dis- 
tribution F(q) peaked at q r with a momentum width 
Sq (x l/l* <C q r and (/>o(p,z) is the ground-state. Taking 
into account (a) and (b), we re- write Eq. (|8| in the form 



:(1-* 2 ) 



dx 2 



2(1 + (3) ) dx 



W(x)F(p)+ 

,2 2/3- 



i 



l + (3 



^(q)(l-X 2 ) 



m 2 LQ 2 
2(1 + 0) J 



W(x)F(p)=0, (10) 



where q = — zV. The lowest energy eigenstates still fulfill 
W(x) ~ 1, but now we must keep explicitly the spatial 
dependence of F(p). 

e 2 F(p) = [£(q) 2 - G(/3)£(q) M (p) + hW z ] F(p) (11) 



^m*o; 2 Vg, which may be calculated as discussed in the 
Letter. 



Density modulations 



Expanding around the roton minimum in the local spec- 
trum, and moving to momentum space (p = zV q ) we 
obtain 



2 F(q) 



. fi 2 (q-q r ) 2 i 2v7 2 



F(q) (12) 



where F(q) is the Fourier transform of F(p). Hence F(q) 
are the eigenstates of the hamiltonian H = - ^ g ~ g ^ 



Note that density modulations are given by the /_ 
amplitudes, in the form: 5n(r,t) = n(r, t) — no(r) cx 
y / no(r)5R(/_(r))., where 5R denotes real part. For the 
uj = case, and since W(x) — 1, then from Eq. ([5| 
one obtains that /_ ~ (E(q)/c)f+. For the uj > 
case, since W(x) ~ 1 still, and we assume that F(q) 
is narrowly peaked around q r , we obtain the same re- 
sult for the trapped case. Hence apart from a constant 
/-(r) cx /+(r) oc ipo(r)F(p). As a result: 5n(r,t) oc 
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no(r)»(F(p)). 

TIME-OF-FLIGHT PICTURES 

Although the most straightforward way of studying 
roton confinement is of course given by in-situ measure- 
ments, time-of-flight(TOF) pictures are expected to show 
a clear difference as well, when compared to the cloverleaf 
pattern characteristic of the phonon collapse (as observed 
in chromium [36] and, more recently, in erbium [3]). More 
specifically, TOF pictures may reveal clear traces of the 
Rashba-like ring characterizing the localized roton-like 
excitations. 

Due to the fast expansion after release, TOF pictures 
are expected to reproduce well the in-situ momentum 
distribution at the time of releasing. We have calculated 
the momentum distribution of the condensate at differ- 
ent stages during the instability dynamics. Our results 
for the same case studied in Fig. 2(a) of the main text 
are depicted in Fig. [5j where we show both the column 
density on the q x -q y plane (i.e. integrated along q z \ and 
that on the q x -q z plane (i.e. integrated along q y ). We 
show in particular the results obtained during the first 
stages of the collapse dynamics. For the parameters of 
Fig. 2(a) of the main text, the Rashba-like ring is clearly 
visible during the first stages of the collapse, being later 
on destroyed due to the subsequent collapse of the spatial 
ring-like density modulations. As a result the ring TOF 
pattern is washed out. In later stages TOF pictures are 
characterized by the appearance of large momentum ex- 
citations on the xy plane (induced by the collapse of the 



local rings) which becomes visible in the q x -q z distribu- 
tion as pronounced jets along x. 

The low-momentum peak always appears saturated in 
our pictures. It was eliminated from Figs. 2b and 2g of 
the main text, since we were only interested in showing 
the Rashba-like feature in momentum space. The large 
low- momentum peak is the result of two effects. The first 
one is purely geometric. The q = peak is always rela- 
tively strong, since it is not spread in the angular variable 
if (as it is the ring feature discussed above). The second 
reason is more physical, and it is linked with the local 
nature of the collapse. Since only a small fraction of the 
atoms is actually participating in the modulational insta- 
bility and in the subsequent collapse, a large fraction of 
atoms remains at low momentum. Note that the low mo- 
mentum peak also broadens due to the generation of low 
energy phonons during the collapse. As a result, the low- 
momentum peak dominates the TOF pictures. This is in 
itself a clear difference with respect to the phonon col- 
lapse, where the initial low- momentum peak transforms 
into a cloverleaf pattern due to the global nature of the 
collapse. 

However, although the relative number of particles pro- 
duced in the Rashba-like ring may be small, the actual 
absolute number of particles may be certainly sizable. In 
Fig. [5] (lower panel) we depict at various times the num- 
ber of particles at different radial momenta. This plot 
shows that a significant number of particles is found in 
the ring before the ring feature is eventually washed-out 
(up to 1500 particles at 21.5 ms). 
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FIG. 5: (Color online) Momentum distribution for different times for the same case depicted in Fig. 2(a) of the main text, 
as it would be revealed by a TOF experiment. In the upper panel we show the q x , q y momentum distribution (integrated over 
q z ) and in the middle panel the q x ,q z momentum distribution (integrated over q y ). The lower panel shows the integrated 
population at different intervals of the radial momentum q — ^/q% + q% • 



